Isotropic-Cholesteric Transition of a Weakly Chiral Elastomer Cylinder 
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When a chiral isotropic elastomer is brought to low temperature cholesteric phase, the nematic 
degree of freedom tends to order and form a helix. Due to the nemato-elastic coupling, this also 
leads to elastic deformation of the polymer network that is locally coaxial with the nematic order. 
However, the helical structure of nematic order is incompatible with the energetically preferred 
elastic deformation. The system is therefore frustrated and appropriate compromise has to be 
achieved between the nematic ordering and the elastic deformation. For a strongly chiral elastomer 
whose pitch is much smaller than the system size, this problem has been studied by Pelcotivs and 
Meyer, as well as by Warner. In this work, we study the isotropic-cholesteric transition in the weak 
chirality limit, where the pitch is comparable or much larger than system size. We compare two 
possible solutions: a helical state as well as a double twist state. We find that the double twist state 
very efficiently minimizes both the elastic free energy and the chiral nematic free energy. On the 
other hand, the pitch of the helical state is strongly affected by the nemato-elastic coupling. As a 
result this state is not efficient in minimizing the chiral nematic free energy. 

PACS numbers: 61.30.-v 61.30.Cz 61.30.Vx 
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I. INTRODUCTION 



An isotropic chiral elastomer can be synthesized 
by crosslinking a chiral nematic polymer melt in the 
isotropic phase. When such a system is brought into the 
low temperature cholesteric phase, the nematic degree of 
freedom orders locally and tends to form a helical struc- 
ture. Due to the nemato-elastic coupling, the polymer 
network tends to stretch along the direction of the lo- 
cal nematic order, which continuously rotates along the 
helical axis. For a system with cylindrical shape, this 
leads to strain deformation which increase linearly with 
the cylinder radius, as illustrated in Fig. 03. Its elastic 
energy cost is formidably high, when the cylinder radius 
is much larger than the helical pitch. This frustration 
due to competition between network elasticity and liquid 
crystalline ordering makes it nontrivial to find the ground 
state of the system in the cholesteric phase. 

This problem was first studied by Pelcovits and Meyer 
[l| using linear elasticity theory. In the limit of infinitely 
strong chirality, it is clear that the system should first 
satisfy the chirality by forming a planar helix along the 
cylinder axis. On the other hand, to avoid large strain 
energy, the solid can only deform uniaxially, which im- 
plies that the nemato-elastic coupling can only be par- 
tially satisfied. Such a state, as illustrated in Fig^2K and 
[2p, is called a planar helix state in reference [l| and a 
transverse cholesteric state in reference [3]. As the chi- 
rality is made weaker, a conical helix state, where the 
precessing director has a nonvanishing component along 
the helical axis, may constitute a better solution. The 
associated solid deformation as well as director pattern 
for this conic state are illustrated in Fig. 03 and 05 re- 
spectively. Warner carried out a nonlinear analysis of 
the same problem using the neo-classical elasticity the- 
ory [3, |3| . Nontrivial dependence of the phase boundary 
on the magnitude of nematic order was identified. The 



multicritical point associated with the planar-conic tran- 
sition (where the first order transition line and the second 
order transition line meet) was also analyzed. 

It is implicitly assumed both in reference [l[ and refer- 
ence m that the pitch of the corresponding nematic liquid 
crystal system (typical < 0.1/xm) is much smaller than 
the system size, e.g. the radius of the cylinder. That 
is, the elastomer is in the strong chirality limit. This 
is certainly correct for many cases. Nevertheless, the 
cholesteric pitch can be continuously tuned by changing 
the concentration of chiral chemical groups during poly- 
merization. In particular it can be tuned to be compa- 
rable to macroscopic length scales, e.g. the system size. 
This is especially true if the system has the shape of a 
thin cylinder or wire. It is therefore interesting and rel- 
evant to study the isotropic-cholesteric transition in the 
weak chirality limit. In this work, we carry out a nonlin- 
ear elasticity analysis of this problem using variational 
methods. We find that in this regime, a double twist 
state has lower free energy than the usual helix director 
pattern. The results obtained by Pelcovits and Meyer 
Warner 0, as well as the authors in this work are 
summarized by the "phase diagram" of a chiral nematic 
cylinder in Fig. [TJ 



II. MODEL 

The total free energy per unit volume of a chiral liquid 
crystalline elastomer crosslinked in the isotropic phase is 
given by 

/ = fel + ,/q, (1) 

where /d is the neo-classical elastic free energy 



(2) 
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FIG. 1: "Phase diagram" of a chiral cylinder in nematic phase. 
The parameters ao, go, and R are defined in Sec.|ITl The tran- 
sition between planar helix state and conical helix state may 
be continuous or discontinuous, while the transition between 
conical helix state and double twist state is expected to be 
discontinuous. 



with 

the deformation gradient matrix defined relative to the 
isotropic reference state f = x, which is subject to the 
incomprcssibility constraint : 

det A = 1. 

As usual, the vector x coincides with the position of the 
mass points in the isotropic reference state and is re- 
ferred to as the Lagrangian coordinate. The vector r 
on the other hand describes the position of mass points 
in the chiral nematic reference state (ground state that 
minimizes the total free energy), and is usually referred 
to as the Eulerian coordinate. As a general property of 
nonlinear elasticity theory, it is important to distinguish 
these two coordinates properly. The symmetric and posi- 
tive definite tensor 1 in the neoclassical elastic free energy 
Eq. ([2]) is called the step length tensor Q of the current 
state, or deformed state @ and describes the statistical 
conformation of polymer chains in the current state. It 
is related to the nematic order parameter Q by 

I = al-6Q. (4) 

where a and b are some microscopic constants. In this 
work we shall always normalize 1 such that it has deter- 
minant one. In the principle coordinate system of the 
nematic order parameter, the step length tensor 1 can be 
represented as a matrix: 

/i \ 

1= i ^iC'-C')nn + C'l, (5) 

V cV 

where ^ is a monotonic increasing function of the mag- 
nitude of the nematic order S, whose detailed functional 



form is irrelevant to our study. For an achiral nematic 
elastomer, C, turns out to be the ratio of spontaneous 
stretch along the direction of the nematic director when 
the system enters the nematic phase from the isotropic 
phase (3]. Due to the incomprcssibility constraint, the 
system shrinks by factor of 1/ in the perpendicular 
directions. Finally we note that in Eq. ^ a constant 
term — 3/Lt/2 is introduced so that the elastic free energy 
vanishes in the isotropic reference state where A = 1 = I. 

The second part /q in Eq. ^ is the Landau-dc Gennes 
free energy for a chiral nematic liquid crystal. Assuming 
that the nematic order is well saturated with fixed mag- 
nitude iS* in the cholesteric state, the relevant nematic 
free energy is the Frank free energy for chiral nematic 
liquid crystals 0, 

/Frank = ^i^i ( V • n)^ + l/sTa (n • V X n - )^ (6) 
-I- ^K3{fi X V X n)^ + K24'V ■ {n-\/h-h\/ -h) , 

where Ki, K2, K^, are splay, twist, and bending constants 
respectively, while q^^ = ^0 is the cholesteric pitch for the 
corresponding chiral nematic liquid crystal. 7^24 is the 
saddle splay constant, whichplays an important role in 
the physics of blue phase [1, U H] . Since the saddle splay 
density is a complete differential, its volume integral can 
be transformed into a surface integral by Gauss' theorem, 
and therefore scales the same as the surface anchoring of 
the nematic director field, which we shall not consider 
in the work. Nevertheless, it is rather straightforward to 
include this surface interaction. Also, it is important to 
note that all the derivatives in Eq. ([6]), Vi = d/dvi, are 
with respect to the Eulerian coordinates, i.e. Cartesian 
coordinates of mass points in their deformed states. This 
is required by the liquid nature of Frank free energy: at 
length scales where the Frank free energy becomes im- 
portant, the system is essentially a liquid. The physical 
quantities of a liquid should be naturally expressed in 
terms Eulerian coordinates, rather than in terms of La- 
grangian coordinates. To avoid confusion in notation, 
we shall use Vi = d/dri for derivative with respect to 
Eulerian coordinates and use da = d/dxa to denote the 
partial derivative with respect to the Lagrangian coordi- 
nates. 

As a first step, let us discuss the total free energy 
Eq. ID) qualitatively. Within the one constant approxi- 
mation of the Frank free energy, and ignoring the surface 
saddle splay term for a moment, there are three natu- 
ral length scales in this problem, q^^ = ^0 is the chi- 
rality pitch, while oq = \/ K/ 11 is the cross-over length 
scale set by the competition between network elasticity 
and nematic director elasticity. The third length scale 
is the radius R of the cylinder. For most liquid crys- 
talline elastomers, we estimate A' ~ 2 - 4 x 10~'^'^N, 
while ^ ~ 10^ — lO^Pa. Therefore ao ~ 1 — lOnm, con- 
stituting the shortest length scale in our problem. On 
the other hand, the chirality pitch (q can vary a lot, typ- 
ically 0.1 fim or smaller for strongly chiral materials but 
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FIG. 2; Elastic deformation and director pattern for various states. A: Reference cylinder in the isotropic phase. B: The 
elastic deformation of the planar helix state at weak chirality limit, studied in Sec. IIVI C: The nematic director pattern of the 
planar helix state. D: The elastic deformation of the conical helix state, studied in references [l|, E: The director pattern of 
the conical helix state, 9 is the conical angle. F: The elastic deformation of the double twist state. G: The director pattern of 
the double twist state. 



may get much larger for weakly chiral materials. In par- 
ticular, it can even be larger than the cylinder radius 
R. Also, the regime a^/R':^ 1 is clearly experimentally 
inaccessible. 

Comparing these three length scales, we are naturally 
lead to the following two distinct regimes: 

1. Weak chirahty regime ao ^ i? ^ ^o- 

2. Strong chirality regime ao,^o ^ R- 

The strong chirality limit has already been analyzed by 
Pelcovits and Meyer as well as by Warner It is 
found that as one tunes the dimcnsionlcss ratio ao/lo = 
aoQo to below a critical value of order of unity, the system 
goes from a planar helix director pattern to a conical 
helix pattern. In this work, wc shall mainly focus on 
the weak chirality limit. Similar to reference [l| and 0, 
we shall use variational methods, proposing two kinds of 
candidate states with certain variational parameters and 
minimizing the total free energy over these parameters. 



III. DOUBLE TWIST OF NEMATIC DIRECTOR 
AND TWIST OF CYLINDER 

Consider a cylindrical block of isotropic chiral elas- 
tomer of radius R, aligning along the z axis. We need 
to find the nematic director field h{x) as well as the clas- 
tic deformation r{x) relative to the isotropic reference 
state that minimizes the total free energy. One possible 
low energy configuration for the chiral Frank free energy 
is a double twist texture, as illustrated in Fig. In 
chiral nematic liquid crystals, the double twist configu- 
ration is energetically favorable if the saddle splay mod- 
ulus is positive and large enough 0,0,11. Accord- 
ing to the current understanding of the blue phase, these 



double twist cylinders pack into a three dimensional peri- 
odic structure with cubic symmetry. In liquid crystalline 
elastomers, due to the nemato-elastic coupling, a double 
twist nematic director texture necessarily induces twist 
of the cylinder, together with a uniaxial stretch A along 
the cylinder axis: 




X, 



(7) 



where 




— sin a z 
cos a z 
1 



(8) 



is a rotation about the z-axis by an angle a z. Using 
the cylindrical coordinate system, the Lagrangian coor- 
dinates of a point x are given by the triplet (p, (jj, z): 



(9) 



In the deformed state Eq. (O, the Eulerian coordinates 
r(x) are given by 




^ cos(0 + az) 
^sin(^-K«z) I , (10) 
A z 



where we have used Eq. ([7]) and Eq. Therefore we 
find 





/ r cos ip ) 




r(a?) = 


r sin (p 











(11) 
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Note that the deformed cyHnder has height L\ and radius 

Let be the unit vector associated with the Eulerian 
cyhndrical coordinate Lp: 

dr 

dip 




(12) 



In terms of the Eulerian coordinates, a double twist tex- 
ture of nematic director is represented: 



e, cosf 



smf 



(13) 



Note that the twist angle e{r) = 0{p/\/X) can be equally 
well represented as a function of Lagrangian coordinate p. 
9{r) satisfies the boundary condition 6{0) — 0, since the 
nematic director is parallel to z on the center axis of the 
cylinder. On the outer surface of the cylinder r = R/\/X, 
9{R/y/X) is free to vary. 

Calculation of the deformation gradient using Eq. ([7]) 
is a trivial and tedious matter. On the other hand, by 
substituting Eq. (fT3|) into Eq. ([5]) we can readily calculate 
the step length tensor 1. Substituting these results into 
Eq. ([2]), we find that the spatially dependent elastic free 
energy density for the proposed double twist solution is 
given by 



/oi — 



1 + A^C^ + aV^C^ 



4C2A 

(C^ - 1) (A^ - a^p 
2a [C^ -\)\^'^psin{2e) 



3C^ + A^ + a^p^ 



(14) 



The spatially dependent Frank free energy density can 
be calculated by substituting Eq. (fT^ into Eq. ([S]), care- 
fully noting that all derivatives are with respect to the 
Eulerian coordinates r. The result is 



Frank 



2 



1 • d0 

— sin 26+ — 
2r dr 

1 • 

— sin 29+ — 

2r dr 



2 r^ 



- K 



24 



sin 29 de 

dr ' 



(15) 



which is identical to that for a chiral nematic liquid crys- 
tal in a double twist cylinder Q. 

In the weak chirality limit, q^R 1, we expect 9{r) 
to be small and linear in r. We can therefore expand the 
elastic free energy density in terms of p [To| and 9{r): 



fci = /o + /2 + higher order terms, 



where 



/o 

/2 



2C 
M 
2C2A 

-2a (C^ 



Cm 3 



{a'p'C' + (C' - 1) 0' (A' - 1) 
1) 9X^/^p) 



(16) 



(17) 



(18) 



are terms of order of r" and respectively. We shall 
ignore all higher order terms in the elastic free energy. 
Note that /o is exactly the free energy density for a 
monodomain nematic elastomer, with anisotropy ratio C, 
undergoing a uniaxial deformation coaxial with the step 
length tensor. Minimizing /o over A we obtain 

A = C ^ /o = 0, (19) 

as expected. Substituting this into Eq. psp . we find 



/2 



2C3 



(C' -1)9- aC 



3/2 



(20) 



which is a complete square. Since /prank is independent 
of a, and since 9 is linear in p as will be shown below, 
Eq. (j20p is minimized by 



a 

f2 



0. 



(21) 
(22) 



Hence there is no elastic free energy cost for the double 
twist state up to the order of {aR)^ . As we shall show 
below, the parameter a is of order of go- Hence (aR)'^ is 
indeed a small parameter in the weak chirality limit. 

Similarly, we expand the Frank free energy in terms of 
9 and r. To the leading order we find 



/prank 

2 \r 



d9 
dr 



K2qo 



9 d9 
r dr 



(23) 

2 9d9 

r dr ' 



which only depends on 9. Note that the bending term 
is of higher order when compared to all other terms that 
we have kept. 

We have to minimize the total Frank free energy den- 
sity 



-Fprank = 27rLA 



/prank rdr 



(24) 



over 9{r) in order to determine the optimal director tex- 
ture. Let us define Rq = R/y/C in order to streamline 
the notation below. Calculating the first variation of 
the Frank free energy, including the boundary terms at 
r = i?^, we find 



SF 



Frank 



^ K2 I dr[ -r9"{r) - 9'{r) + 



'1] 



S9{r) 



2ttLX " Jq \ ^ ' ^ ' r J 
+ [K2{ 9{R^) + Rc9'{R^) - R^qo ) - 2K2i9{R^)] 69{R^). 

(25) 

Since the twist angle 9{r) is free to vary on the bound- 
ary r = Ri^, we have to set both the integrand and the 
boundary term to zero in order to find the minimizing 
solution. This leads to the following two Euler-Lagrange 
equations 



r9"{r) 



r{r) + ^ = 0, 
r 



(26) 



i9{Rc)+Rc9'{Rc)-Rcqo)~2rj9iR^)=0, (27) 
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where r] = K24/K2 is a dimcnsionlcss ratio. Solving these 
two equations we find 



0ir) 



qo 



qop 



(28) 



which exphcitly shows that 6{r) is indeed Hncar in r. The 
twist angle on the boundary is given by 



d{Rc) 



2(l-r;)yC' 



which serves as a small parameter controlling the validity 
of the pcrturbativc analysis. Substituting Eq. (pS)) into 
Eq. (1211) we find the parameter a given by 



a 



(C^-l)go 

2(1 



(29) 



which is indeed a constant, of the same order of go, and 
independent of r. Substituting Eq. ((28|) into Eq. ((23)) we 
find the Frank free energy density, which is also the total 
free energy density (since the elastic free energy vanishes 
at the order of (ai?)^), to be given by 



/tot — /] 



Frank 



K2ml 

'2(l-r,)- 



(30) 



Summarizing Eq. Eq. ^ and Eq. we find 

that if ?7 < 1, our perturbative calculation is quantita- 
tively good in the weak chirality regime where 90^/2(1 — 
??)-\/C <C 1. The double twist state is very efficient in min- 
imizing both the elastic free energy and the Frank free 
energy. In particular, when the saddle splay constant 
K24 vanishes, 77 = 0, and therefore the total free energy 
Eq. ([50)) also vanishes. Note that the total free energy is 
positive definite if 7/ = 0. Hence the double twist state 
is clearly the ground state, at least up to the order of 
(q^Kf'. By contrast, for a cholcsteric liquid crystal with 
K2A = 0, the blue phase is clearly not the lowest energy 
state, compared to the usual helical state. This shows 
that unlike in the blue phase of cholesteric liquid crystal, 
the saddle splay constant if 24 does not play an impor- 
tant role in the formation of the double twist pattern in 
a cholesteric elastomer. When qoR/2{l — ?/)\/^ is compa- 
rable or larger than unity, the higher order terms of the 
elastic free energy and the Frank free energy can not be 
neglected, and one has to minimize the full free energy 
Eq. (fT4)) and Eq. (fT5)) . Finally if 77 > 1, a perturbative 
calculation in power of aR is qualitatively incorrect, no 
matter how small the parameter qoR is. We must mini- 
mize the full elastic free energy Eq. (fH)) and Eq. (fT5)) . 



IV. HELICAL STATE 

In the weak chirality regime that we are interested in, 
<Zo«o QoR ^ 1: the elastic energy scale (per unit vol- 
ume) fi is much larger than the chiral Frank energy scale 



Kq^. Therefore the conical helix state studied in refer- 
ence [l| and can never be the ground state, as it only 
partially minimizes both the Frank free energy and the 
elastic free energy. There is however, another potential 
candidate for the ground state, which can minimize the 
elastic free energy up to the leading order. Let us con- 
sider a planar helix director pattern along the cylinder 
axis, where the nematic director remains perpendicular 
to the cylinder z-axis and rotates around this axis with 
pitch a 



cylinc 
lul: 

h{z) 



cos a z + iy sin a z = (a z) e^; 



(31) 



In the following, we shall use both dyadic notation and 
matrix notation of tensor quantities. The corresponding 
local step length tensor is given by 



where 



l(z) = 0^[az)\[z = 0)O3(-az), 



l(z = 0) = 




(32) 



(33) 



is the step length tensor at the plane 2; = 0. This vari- 
ational form of nematic director field is the same as the 
planar helix sate considered in reference [l| . 

Due to the nemato-elastic coupling, the polymer net- 
work prefers to stretch along the local nematic director. 
This however implies that the direction of local strain 
deformation rotates by an angle 7r/2 between two cross 
sections ^o/4 apart along the cylinder. This leads to an 
additional strain energy density /.t (ai?)^ that is quadratic 
in the cylinder radius. For a fat cylinder (or in the 
strong chirality limit), aR ^ 1 and this strain energy 
is prohibitively high |12| . For a thin cylinder (or in the 
weak chirality limit), however, aR <^ 1 and this addi- 
tional strain energy only constitutes a perturbation to 
the strain energy of the corresponding uniform deforma- 
tion. Nevertheless, to reduce the additional strain energy 
at the order of fi (aR)'^, the system can globally twist in 
the direction opposite to the nematic helix. The over- 
all nonuniform deformation, shown in Fig. 03, is repre- 
sented by the Eulerian coordinates as functions of the 
Lagrangian coordinates: 



f{x) 



= Oz{az) 



'■)■ 


X 




A 











1 











1 

Vx 



(34) 



0,(-(a + /3)z) 



where (3 measures the global twist of the solid. We also 
note that due to the additional inhomogeneous strain de- 
formation, the inverse pitch of the director helix a is 
generically different from the inverse pitch go of the cor- 
responding liquid crystal system. This will be made clear 
in the calculation below. 

It is important to note that the matrix A defined 
through Eq. ([M)) explicitly depends on the coordinate 
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z and is not the deformation gradient A^a- The lat- 
ter should be obtained by taking partial derivative of 
Eq. ([34|) with respect to Lagrangian coordinates x,y,z. 
Being derived in this way, the deformation gradient ma- 
trix naturally satisfies the following compatibility condi- 
tions: 

Substituting the deformation gradient and nematic or- 
der parameter Eq. into Eq. integrating over the 
reference volume of the cylinder, and dividing it by the 
total volume ttR^L, we obtain the clastic free energy den- 
sity for the proposed deformation gradient as 



/el [a,/?, A] 



/o 

/2 



/O + /2, 



" A 

.2 (C^- 



(C^ + A^) (3^-2 (C^ 



1) (a3/2 - 1 



(35) 
(36) 

(37) 



Note that /□ is independent of the cylinder radius R and 
is identical to the free energy of an achiral nematic elas- 
tomer undergoing uniaxial and homogeneous deforma- 
tion. By contrast, /2 is proportional to R^, and quadratic 
in a and /3. /2 is clearly due to the inhomogeneous de- 
formation. The ratio between /o and /2 scales as (aR)^ 
as discussed earlier. The dimensionless ratio aR char- 
acterizes the importance of chirality in this problem. In 
reference i and i this ratio is implicitly taken to be 
large at the very beginning. In this work, we shall as- 
sume it to be a small number. More precisely we shall 
assume another dimensionless ratio qo R <SC 1. Also we 
shall see below that for this proposed variational solu- 
tion, the Frank free energy scales the same as /2, hence 
it is reasonable to first minimize /q and then the sum of 
/2 and /prank- Minimization of /o leads to 



A 

/o 



0. 



(38) 
(39) 



That is, the local elastic deformation is identical to that 
of a homogeneous achiral nematic elastomer. Inclusion 
of /2 and /prank iuduccs small change of A at order of aR. 

The Frank free energy density for the proposed director 
pattern Eq. ([?T|) can also be easily calculated. Again we 
have to be careful with the derivatives in Eq. ([6]) that are 
with respect to the Eulerian coordinate r. After some 
tedious calculation and replacing A with wc find 



Frank 



i A-2a2C - K2aqo VI + \K2ql , (40) 



which is independent of elastic constants i^i, and 

We still need to minimize the sum of /2, given in 
Eq. psp and /prank Eq. (j40p , over the remaining two vari- 
ational parameters a, (3. Since /prank does not depend on 



/?, we minimizes /2 over (3 and find 



f2 



2(3/2 
R^a^ (C'-l)V 



16(3 



(41) 
(42) 



We note that as long as ( 7^ 1, and a 7^ 0, the global 
spontaneous twist of solid /3 does not vanish. More im- 
portantly, unlike the double twist state, the planar helix 
state considered here does cost elastic free energy at the 
order of (aR)^. 

We can now minimize the sum of Eq. and Eq. 
over a, which leads to 



i?2 (^3 „ lY 



(43) 



Remembering uq/R < 10~^ even for R = 1mm, and 
C 7^ 1, the first term in the denominator can be safely 
ignored and a can be approximated as 



i?2 ((3 



1) 



ooV 

r) 



qo < qo- 



(44) 



This result indicates that the helix of nematic director is 
strongly resisted by the nemato-elastic coupling energy 
and the pitch becomes much longer than the correspond- 
ing value £0 in the nematic liquid crystals. Substituting 
Eq. (glD into Eq. ^ and Eq. we find 

The total free energy is therefore positive, in strong con- 
trast with the double twist state we considered in the 
preceding section. The planar helix state considered here 
is therefore not efficient in energy minimization. This is 
clearly due to the extra elastic free energy cost Eq. 
caused by the nemato-elastic coupling. 



V. DISCUSSION AND CONCLUSION 

We have shown in this work that in the weak chirality 
limit qQOQ <ti qoR ^ 1, the double twist state minimizes 
both the Frank free energy and the elastic free energy 
up to the order of (qoR)'^, and is therefore a good can- 
didate for the real ground state. The planar helix state, 
on the other hand, is strongly influenced by the nemato- 
elastic coupling, with the pitch much longer than the cor- 
responding value in cholesteric liquid crystal. As the di- 
mensionless parameter q^R becomes larger than one, the 
elastic energy cost due to inhomogeneous strain, scaling 
as iJ,{qoR)^, dominates all other terms. When qoR ^ 1 
and qoOo ^ 1, the ground state is likely to be the conical 
helix state with 9 « 7r/2, according to th e studies in 
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references [l|, [2] ■ The conical state and the double twist 
state are qualitatively different, and can not be mutually 
accessed in a continuous fashion. Therefore the afore- 
mentioned two regimes must be separated by a first or- 
der phase transition, located around qoR ~ 1. Study of 



this transition is technically challenging and is beyond 
the scope of this work. 
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